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ABSTRACT 

We present the results of a Spitzer photometric investigation into the IR 

excesses of close binary systems. In a sample of 10 objects, excesses in IRAC 

and MIPS24 bands implying the presence of warm dust are found for 3. For 2 

objects we do not find excesses reported in earlier IRAS studies. We discuss the 

results in the context of the scenario suggested by Rhee and co-workers, in which 

warm dust is continuously created by destructive collisions between planetary 
O 

in 
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bodies. A simple numerical model for the steady-state distribution of dust in one 
IR excess system shows a central clearing of radius 0.22 AU caused by dynamical 
perturbations from the binary star. This is consistent with the size of the central 
clearing derived from the Spitzer spectral energy distribution. We conclude that 
close binaries could be efficient "destroyers of worlds", lead to destabilize the 
orbits of their planetary progeny by magnetically-driven angular momentum loss 
and secular shrinkage of the binary separation. 

Subject headings: binaries: close - - circumstellar matter - - infrared: stars - 
stars: individual (AR Psc, II Peg, UX Ari) -- planets and satellites: dynamical 
evolution and stability 



1. Introduction 

Perhaps the most spectacular discovery of the Infrared Astronomical Satellite (IRAS) 
was the detection around a small handful of normal main-sequence stars of excess infrared 
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(IR) radiation from remnant disks now understood to represent debris related to planet 



formation (e.g. lAumann et al.lll984[ ). These disks are largely dissipated in in a few hundred 



Myr by collisional degradation of larger bodies combin ed with Poynting- Rob ertson (PR) and 
stellar wind drag, radiation pressure and evaporation (lAumann et al.lll984f ). Recent Spitzer 
surve ys at 70 /xm find a disk fraction of ~ 10% fo r F, G and K stars with ages up to 3 Gyr 
or so farvden et al.lbood : Eillenbrand et al-lbood ). 



One other group of stars with IR excesses has remained inconspicuous. The RS CVn 
class of tidally-interacting close binaries, in which at least one component has generally 
reached the sub-giant evolutionary phase (JHalllll976l ). were discovered by IRAS to have an 
IR excess occurren ce rate of about 40% — twice as high as main -sequence stars and 3-4 times 



higher than giants (IBusso et al.lll987l . ll988l ; IScartriti et al.lll993l and references therein). The 
observed excesses typically amounted to several percent of the stellar bolometric luminosity 
and corresponded to a blackbody temperature of order 1500 K. Within current scenarios 
explaining IR excesses around late type stars, RS CVns should not have any: th ey represen t 
a significantly older population than the few 100 Myr dusty disk lifetime (e.g. lEkerl 119921 ). 
and have not undergone any evolutionary phase during which substantial mass loss and 
dust formation is expected. Despite this puzzle, the IR excess RS CVns have received scant 
attention since their initial discovery. 



star" (jZuckerman et al. 



More recently. IWeinbergerl (120081 ) found BD +20 307, "the dustiest known main-sequence 



20081 ) . to be a close spectroscopic binary with a period of 3.4 days 
and estimated its age to be greater than 1 Gyr. The parallel between the 3.4 day binary 
BD +20 307 and the IR excess RS CVns — spectroscopic binaries with periods of a few 
days — is striking. The dust emission of BD +20 307 amounts to an unusu ally large 4% of it s 
bolometric l umin o sity, and has an un usual ly high temperatu r e of ~ 650 K (jSong et al.ll2005l ). 



Song et al.l ( 120051 ) , iRhee et al.l ( 120081 ) , and IZuckerman et al.l ( 120081 ) suggest the warm dust in 
BD +20 307 likely arises from recent pulverising collisions between planets or planetesimals. 
That such close binaries can harb or planetary systems has been demonstrated by the planet 
around HW Vir (JLee et al-lbooih . 



Here we present a Spitzer photometric survey of a small sample of 10 RS CVns to probe 
deeper into the purported IR excesses of this class. 



2. Observations and Data Analysis 



The data presented in this paper were obtained by the Spitzer Space Telescope (IWerner et al. 



20041 ; iGehrz et al.l 120071 ) as part of the General Observing program PID 20713. Infrared Ar- 
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ray Camera (IRAC; Fazio et al.ll2004i ) observations were made in High Dynamic Range mode 
between 2005 Jul and Dec. Targets were observed in 3.6, 4.5, 5.8 and 8.0 /zm bands for 12s 
in each position of a 12 point Reuleaux pattern dither map, for a total e ffective integration 
time of 124.8s. Multiband Imaging Photometer (MIPS; iRieke et al.ll2004l ) observations were 
obtained in the 24, 70 and 160 /xm bands in photometry mode between 2005 Nov and 2007 
Jan, with an exposure time of 10s, repeated for 2 cycles (24 /xm) or 3 cycles (7 and 160 /jm). 



Redu ced image mosaics were produced with the IRACproc custom software (ISchuster et al. 



20061 ). while for MIPS we used the Spitzer Science Center pipeline mosaics. These were 
analysed using standard aperture (MIPS) and Point spread Function (PSF) base d (IRAC) 
photometry, the latter adopting the procedure described in iMarengo et al.l (120091 ). that en- 
able deriving the flux of core-saturated sources. All photometry results are reported in 
Table [U 

Only one source (UX Ari) was detected at 70 [im at > 3a. For two other sources (CF 
Tuc and II Peg) we found 2<r level signals. For all other sources we report 3er upper limits. 
The well known MIPS 160 micron near-IR spectral leak prevented us from using data taken 
in this band. 



Our analysis procedure was similar to that followed by lBusso et al.l (ll987U1988l ); IScartriti et al 



(119931 ) and involved comparison of stellar model fluxes, fitted to observed optical and near- 
IR photometry, with the Spitzer IR photometry. For each binary listed in Table (H we 
modelled BVRIJHK photometry culled from literature with a spectral energy distribution 
(SED) model comprising the sum of two Kurucz photospheric modelsj (one for each stel- 
lar component in the case of double-lined spectroscopic binaries, weighted by stellar radius 
and effective temperature according to r 2 T e 4 ff ). Since surface spots resulting from magnetic 
activity on RS CVns can give rise to detectable secular photometric changes (thus leading 
to some having been designated as variables), we adopted photometric uncertainties that 
generally encompassed observed ranges of variation s for e ach s tar. For dwarfs stars, we 
adopted the effective temperature-radius relation of iGrayl (119921 ) with initial values for ef- 
fective temperature from the literature; for more evolved components we culled values for 
temperature, radius and/or mass and surface gravity from the literature. Final radius and 
effective temperature model parameters for each system were honed by exploring the grid of 
Kurucz models using a \ 2 minimisation. 

Final stellar parameters are listed in Table |2j Comparison of adopted Kurucz models 
with Spitzer photometry revealed satisfactory fits for all but the three sources: AR Psc, II Peg 
and UX Ari. For these systems, no combination of Kurucz photospheres could reproduce the 



1 http://wwwuser.oat.ts.astro.it/castelli/grids.html 
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measured IRAC and MIPS fluxes: in all three cases, the observed IR fluxes were too high 
compared to the predicted model fluxes. Including a blackbody as an additional component, 
fitting the temperature and normalization as free parameters (see Table [2]), greatly improved 
the SED model fits. Results are illustrated in Figure [TJ, together with IRAS Fain t Sour ce 
Catalog v2.0 12 \xm fluxesj. In agreement with earlier studies (e.g., IScaltriti et al.l ( 119931 )). 



the best-fit blackbody temperatures lie in the range 1300-1900 K with an uncertainty of 200- 
300 K. The low-level 70 \im emission in II Peg is consistent with the long wavelength tail of 
the hot blackbody excess, while that for UX ARI hints at the possible presence of additional 
cold dust. For other binaries that do not show a detectable IR excess, we have estimated 3<r 
upper limits. Assuming a nominal blackbody temperature of 1500 K, we find that for these 
systems the blackbody luminosity must be less than « 2 — 2.5% of the bolometric luminosity 
(see Tabled]). 

To estimate more realistic dust temperatures, non-blackbody circumstellar dust with 
an emissivity law of the type e = 1, A < Ao; e = (Ao/A) 9 , A > Ao, where Ao corresponds 
approximately to the grain size and q is a parame ter of order unity, depending on grain type, 



as described by, e.g., iBackman fc Parescd (119931 ) were used. In general, fits using a grain 



emissivity law (Table [2]) required lower temperatures. For both II Peg and AR Psc, the 
best-fit values of Ao and q are 3 jum and 1 respectively, whereas for UX Ari these parameters 
were not well-constrained by the data. 

We note two systems in Table [T] for which we find no significant excess at levels below 
earlier IRAS-based detections; either the earlier detections were spurious, or the excesses 
have waned in the intervening years. 



3. Discussion 

The span from 3.6-24 /mi covered by our Spitzer IRAC and MIPS photometry is very 
effective for diagnosing the presence of warm and even hot dust. The requirement for a 
purely stellar SED is that the IR photometric indices lie along the Rayleigh- Jeans tail of 
the photospheric IR (essentially blackbody) continuum. The SEDs illustrated in Figure [1] 
fail to do this in the sense that the observed slope is incorrect: the departure from stellar 
SEDs increases toward longer wavelengths. From the blackbody fits, we have determined 
the fraction of the stellar bolometric luminosities reprocessed into IR emission and list these 
in Table[T] we find values ranging from 0.7 to 1.9%. 



2 http://irsa. ipac.caltech.edu 
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Our det ection of IR excess in the spectrum of II Peg c onfirms the earlier photometric 
detection by lBusso et a l. (1987 , see also IScaltriti et al.lll993l) and the preliminary ISOPHOT 
spectrum obtained by lRodono et al.l (119981 ) . Of our other two detections, UX Ari w as found 
to have a substantial near-IR excess in JHK photometry by 
not found to have a significant excess in IRAS photometry by 
instead represents a new IR excess detection. 



Verma et al.l (119831) but was 
Busso et all Jl98sh . AR Psc 



As noted by IScaltriti et al.l (119931 ) , the observed IR excesses cannot be fitted with a 
standard free- free emission spectrum from, e.g., an ionized wind, whose u os frequency de- 
pendence (e.g. lWright fc Barlowlll975l ) is much too flat. Moreover, the observed fluxes would 
imply mass l oss rates of o r der 10 ~ 8 M& yr -1 — orders of magnitud e higher than expected for 



these stars. iBusso et al.l (119871 . Il988l ) and IScaltriti et al.l (119931 ) interpreted the observed 
excesses in terms of dust shells. While the implied dust temperatures are high in the context 
of dust in debris disks, we emphasise that unresolved low-mass stars or brown dwarfs cannot 
explain the excesses be cause typical lumin osities of these objects for the required tempera- 
tures are L < 10~ 3 L & (jRiddick et al.ll2007l ) which is more than an order of magnitude lower 
than observed. 



Busso et al.l (119871 ) and IScaltriti et al.l (119931 ) suggested the dust arises from the accu- 
mulation of stellar winds driven by magnetic activity. However, in order to reach tempera- 
tures of ~ 1200 — 1400 K, dust must reside close to the central stars at typical distances of 
~ 0.1 — 0.3 AU . Dust in such a loca t ion is short-lived and any grains that form will be rapidly 
dissipated by PR drag. iRhee et al.l (120081 ) estimate the lifetimes of larger grains of a few //m 
in size around BD +20 307 to be only a few decades. Minimum dust masses able to explain 
the excesses of II Peg, AR Psc and UX Ari are 1.3 x 10" 12 , 2.8 x 10~ 13 and 4.7 x 10~ 13 Mq. 



respec tively; a standard gas:dust ratio implies total masses 100 times higher. IWood et al. 
(120051 ) found very active stars might lose mass at rates perhaps up to 1O~ 12 M yr _1 (100 



times that of the Sun), meaning 10-100 years of accumulated wind would be necessary to 
explain the dust. It seems highly unlikely that such fully-ionized, magnetically-driven winds 
and coronal mass ejections with outflow velocities of several hundred km s _1 would accumu- 
late in this way. 



We are lead to the suggestion of lRhee et al.l (120081 . see also lWeinbergerll2008l : IZuckerman et al. 



20081 ) that such warm dust originates from collisions between planetary bodies. Since grains 
produced in this way will be dissipated relatively quickly, such collision events must be 
recent. 

The problem, then, is why stars of such ages as the RS CVn-type binaries should be 
producing dust when debris disks around single star s and wide binaries are la rgely dissipated 
in a few hundred Myr. In the case of BD+20 307, IZuckerman et al.l (120081 ) noted possible 
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systemic radial velocity variations and suggested an unseen brown dwarf or low-mass star 
could have destabilized the orbits of rocky bodies. Such destabilization, occurring many 
millions of years after a "first generation" debris disk would have been mostly dissipated, 
could stimulate new dust production in Gyr old systems, reviving whatever debris disk 
remains around these stars. 

We suggest another possible mechanism resulting from the conspicuous property of 
tidally-interacting binaries not shared by their non-interacting single and binary brethren 
of similar age: copious magnetic activity manifest in the form of high levels of surface 
spots, chromospheric and X-ray emission. Such activity is thought responsible for the rota- 
tion al braking thr ough angular momentum loss of main-sequence stars from zero-age onwards 
(e.g. lMestellll984[ ). Angular momentum loss of the dwarf precursors to RS CVns, such as 



BD+20 307, as well as the more X-ray bright RS CVns themselves, results in substantial 
orbital shrinkage (e.g. by factors of up to 2 or more) through spin-orbit tidal coup ling; on 
timescales of hundreds of Myr to a Gyr (e.g. lEggleton fc Kiseleva-Eggletonl |2002| ). This 
change in the orbit causes proportional changes in the locations of both mean motion reso- 
nances and secular resonances in whatever planetary system the binary has. We postulate 
that the sweeping of these resonances can lead to both destabilization of the orbits of small 
bodies and resonant capture; both of these processes can increase the collision rate among 
these bodies. This mechanism would provide a more general means of stimulating new eras 
of copious dust production and debris disk revival in Gyr old systems without the need to 
invoke perturbations from additional companions. 

One potential problem for a collisions! o rigin of the dust is the h igh frequency of 
detection — 40% according to lBusso et al.l (119881 ) and IScaltriti et al.l (119931 ). Our pilot study 
sample is partially biased toward systems with previously identified excesses. If our non- 
detection of excesses in CF Tuc and WY Cnc imply earlier detections were spurious, the 
excess frequency implied by this and our detection of one excess in five objects not previ- 
ously studied is ~ 20%. An evolutionary timescale of 1 Gyr then implies an integrated warm 
dust production period — whether in a single or multiple episodes — of 200 Myr. For a 100 yr 
dissipation timescale of 10~ 12 M of dust, the total mass required is ~ 2 x 10~ 6 M Q or about 
1 Earth mass, consistent with a planetary origin. 



3.1. Dynamical Debris Disk Model 

To examine the inter action between a close binary star and a debris disk we ran a 
numerical model (see, e.g. IWilner et al.ll2003l ) of a dust cloud around II Peg using the dust 



properties consistent from our Spitzer photometry. II Peg was chosen because it is the best 
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studied of the three systems for which we detect IR excesses: we assumed masses of 1.67 and 
O.23M , a semi-major axis of 0.087 AU, ec centricity of 0.0, and t otal l uminosity of 4.87L ffi 
based on the stellar parameters reported by lOttmann et al.l (119981 ) and lFrasca et al.l ( 120081 ). 



For this initial model, we focus on dynamic al clearing and neglect g rain-grain collisions, 
which can also alter the disk morphology (e.g. IStark fc Kuchnerl 120081 ) . Dust was released 
from 1 AU with initial eccentricities and inclinations uniformly distributed over the ranges 
0-0.1 and 0-0.05, respectively, with the rest of the orbital elements uniformly distributed 
over 0-27T. We numerically integrated the orbits of 620 dust grains, subject to gravity from 
the stars, radiation pressure, and PR drag, until the particles either collided with one of the 
stars or was ejected beyond 4 AU. Another parameter in these models is 0, the radiation 
pressure force on the grain divided by the gravitational force. In the vicinity of a star with 
the luminosity of II Peg, a 3 /iin radius silicate sphere consistent with our models for the 
photometry corresponds to j5 = 0.463. We chose this value, and neglected the slight radiation 
pressure from the faint secondary star; we leave this detail for future models. At regular 
time intervals, we accu mulated the particle positions in a histogra m, in the manner of many 
previous authors (e.g. IWilner et al.l 120031 ; IStark fc Kuchnerl 120081 ). This histogram models 
the steady-state surface density of the cloud, which is illustrated in Figure 

Figure [2] shows a disk with a central clearing, centered on the barycenter of the binary. 
White dots indicate the locations of the two stars. The dust surface density distribution 
i s azim uthally symmetric, consistent with the analytic predictions of iKuchner fc Holman 
(120031 ). There are 2-3 ring-like density enhancements, which might be artifacts of the narrow 
particle size distribution in these initial models, but no signs of clumps or spiral features. 
For interpreting the Spitzer data, the most important result is the size of the central clearing. 
The radius of this clearing in the dynamical model is 0.22 AU and similar to the radius at 
which first order mean motion resonances overlap (and actually coincident with the 4:1 mean 
motion resonance). This is consistent with our models for the spectral energy distribution for 
which the grain emissivity law temperature implies a radius 0.241q 06 AU. A dust sublimation 
temperature of 1700 K implies that the dust can survive at 0.22 AU, and that the clearing 
radius is determined by the system dynamics, rather than grain evaporation. 

In Section &\ we noted that IRAS excesses for two systems i n our Spitzer sample are 
not confirmed. The Poynting- Robertson time for 3 /im grains is (jWyatt fc Whipple! Il950[ ) 
T PR = 400(a /lAU) 2 L Q /f3L* yr. For the grains in the II Peg model, a = 1AU, L* = 
4.87L Q ,/3 = 0.463, and the PR time is 177 years. But if the grains are launched near the 
star at 0.3 AU instead of 1 AU, the PR time drops to 16 years — sufficiently short for IRAS 
epoch dust to have dissipated by the time Spitzer observations were made. Either the earlier 
detections were spurious, or excesses could have waned, suggestive of dust being produced 



in a succession of discrete events. 



4. Conclusions 

A Spitzer survey of ten tidally-interacting binary systems finds IR excesses for three. 
The excesses are consistent with significant amounts of warm dust close to the central stars. 
A numerical model for the steady-state distribution of dust in II Peg finds a clearing radius 
of 0.22 AU, consistent with the observed optical-IR SED, and suggests the warm dust tem- 
peratures result from inward migration of dust due to PR drag and central clearing caused 
by dynamical perturbations. 

Reported IRAS excesses for CF Tuc and WY Cnc are not confirmed by Spitzer, indicat- 
ing a possible waning of dust in the last few decades. This is consistent with PR drag, for 
which we find the timescale for dust depletion around II Peg could b e as short as a deca de or 



so, suggesting ongoing dust production. We echo the suggestion of iRhee et al.l (120081 ) that 
dust could arise from collisions of planetary bodies, and speculate that these old binaries 
can enter epochs in which they destabilize the orbits of their planets through shrinkage of 
their own binary orbit due to magnetically-driven angular momentum loss. 



M. Matranga was supported by Spitzer Contract 1279130; JJD and VLK were funded 
by NASA contract NAS8-39073 to the Chandra X-ray Center. Observations were made with 
Spitzer, operated by JPL under NASA contract 1407. 
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Table 1. Summary of photometric data for the program stars 



Target 



3.6 /im 



4.5 /im 



MAGNITUDE 
5.8 fim 8.0 fim 



24 fim 



70 fim 



EXCESS 
t /Uol [lO- 2 ^ Prev. work 



AR Lac 15 - 3 
AR Psc 1 ' 2 
CF Tuc 3 - 4 
II Peg 2 ' 5 - 6 
TY Pyx 7 
UV Psc 8 
UX Ari 2 > 3 - 9 
V471 Tau 10 > n 
WY Cnc 12 - 13 
XY Uma 14 



4.23±0.02 4.28±0.03 4.24±0.03 4.25±0.03 4.07±0.08 > 3.06 

4.82±0.03 4.87±0.03 4.83±0.03 4.81±0.03 4.69±0.05 > 3.26 

5.38±0.03 5.42±0.03 5.41±0.04 5.37±0.04 5.20±0.05 4.65±0.51 

4.53±0.03 4.56±0.03 4.54±0.04 4.51±0.04 4.28±0.04 3.75±0.56 

5.20±0.03 5.27±0.03 5.26±0.03 5.26±0.04 

7.22±0.04 7.19±0.04 7.22±0.05 7.13±0.05 6.99±0.07 > 4.50 

3.92±0.03 3.91±0.02 3.87±0.02 3.89±0.02 3.54±0.04 3.09±0.30 

7.19±0.03 7.20±0.03 7.18±0.04 7.15±0.04 7.05±0.07 > 2.55 

7.44±0.03 7.45±0.04 7.41±0.05 7.40±0.05 7.31±0.07 > 2.97 

7.16±0.02 7.17±0.03 7.15±0.05 7.11±0.04 7.13±0.07 > 3.42 



< 2.0 

1 5+ ' 4 
1, °— 0.4 

< 2.3 

i q+0.5 

<1.7 

< 1.9 

7 +0 ' 2 
u -'-0.2 

< 2.4 

< 2.2 

< 1.9 



Y" a, 6 

N d 

Y e - N E 
Y b 



INote. MIPS magnitudes from aperture photometry using the following source radii and sky background annuli. 70(im: 

16" source, 18-39" background, aperture correction factor 2.04. 24/i.m: 13" source (6" for V471 Tau), 20-32" background, aperture 
correcti on factor 1.17 (1.70 for V4 71 Tau). — 'Me asu red dust luminosity v alues and 3a upp er l imits — B VRIJHK p hot ometry refer- 
ences: l ICutispoto et alj feOOlr i ; 2 lMessinal [1200811 ; 3 lBusso etal] lll988f); 4 ICutispotol 11199111; 5 IScaltriti et"ai] Jl99 j) ; 6 iBvrne et al 



1 1996 ); 7 lBusso et ~~dl987h; " 8 lKiurkchieva et alj <2005h ; 9 lAIekseev fc Kozhevnikoval | |2004.1; 10 
l 2005h ; 12 iMilano et al.1 jl986h: 13 JHeckert et alj Jl99iHi; 14 lKiurkchieya et alj feOOCft 15 iRodono e i al 



h983h 



1 1 .amo elu^ 19891) ; " llbanoglu eHd 
. (1986) — Previous work ex- 



, |1V111CX1H_> LL CU. | -L^7 0UIJ. MCUV^lh LI Cll. (±^31?/, XV | Lll JVL111U Vd UL dll. \Z,VJUUy UUUUI IU Gb CL1.| U_L^7UU|f X 1LV1UU3 WU1A GA~ 

from the literature: HBusso et al.1 lll988t) ; HScaltriti et al.Nl993h ; c lRodono et alJll99Sl) ; d lBusso et al.1 dl987h ; e IVerma et"ail 
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Table 2. Stellar parameters and best-fit model results 



BEST-FIT MODELS 



Target 



Photospheres 
Sp. Type a P [day] a T eff [K] b R 2 /Ri 



[K] 



Dust 
T 



[K]d 



AR Lac 


G2IV/K0IV 


1.98 


5750/4750 


2.4 






AR Psc 


K1IV/G7V 


12.25 


4500/5250 


0.8 


1900+ 300 
iyuu_ 300 


UOOJjgj 


CF Tuc 


G0V/K4V 


2.80 


6500/4750 


3.1 






II Peg 


K2IV/M2V 


6.71 


4500 




1700±|gg 


1200+ 200 


TY Pyx 


G5IV/G5IV 


3.20 


5750/5500 


0.6 






UV Psc 


G5V/K2V 


0.86 


5500/5250 


1.1 






UX Ari 


G5V/K0IV 


6.44 


5250/4750 


3.6 


1300+ 200 




V471 Tau 


K2V/DA 


0.52 


5500 








WY Cnc 


G5V/K0V 


0.83 


5500/3500 


0.8 






XY Uma 


G5V/K5V 


0.48 


5750/4250 


0.7 







Note. -- a Spectral types and periods from IStrassmeier et al.1 (119931 ) 

- b Based on Kurucz models; the error associated with the best-fit 
photospheric temperatures is ~ ±250K - c Blackbody dust temperature 

- d Grain emissivity law dust temperature (see text). 
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Fig. 1. — Observed and model SEDs for the IR excess stars found in this survey. The black (dashed) line 
represents the photospheric component, the red (dot-dashed) line corresponds to the blackbody dust compo- 
nent and the blue (solid) line represents the sum of photospheric and blackbody dust models. Photometric 
measurements are represented by different coloured symbols: green - BVRIJHK, red - IRAC (II, 12, 13, 
14) and MIPS (M24, M70). IRAS 12 fim data are indicated by hollow black circles. 
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Fig. 2. — The steady-state surface density distribution of a dust cloud model for II Peg. See 
Section 3 for details. White dots indicate the locations of the two stars. The dust surface 
density distribution is azimuthally symmetric and shows a central clearing centered on the 
barycenter of the binary. The radius of this clearing in the dynamical model is 0.22 AU, which 
is consistent with grain temperatures derived from the observed spectral energy distribution. 



